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ABSTRACT 

We investigate evolution of the intracluster medium (ICM), considering the 
relaxation process between the ions and electrons. According to the standard 
scenario of structure formation, ICM is heated by the shock in the accretion flow 
to the gravitational potential well of the dark halo. The shock primarily heats 
the ions because the kinetic energy of an ion entering the shock is larger than 
that of an electron by the ratio of masses. Then the electrons and ions exchange 
the energy through coulomb collisions and reach the equilibrium. From simple 
order estimation we find that the region where the electron temperature is 
considerably lower than the ion temperature spreads out on a Mpc scale. 

We then calculate the ion and electron temperature profiles by combining 
the adiabatic model of two-temperature plasma by Fox & Loeb (1997) with 
spherically symmetric N-body and hydrodynamic simulations based on three 
different cosmological models. It is found that the electron temperature is 
about a half of the mean temperature at radii ~ 1 Mpc. This could lead to 
an about 50 % underestimation in the total mass contained within ~ 1 Mpc 
when the electron temperature profiles are used. The polytropic indices of the 
electron temperature profiles are ~ 1.5 whereas those of mean temperature 
~ 1.3 for r > 1 Mpc. This result is consistent both with the X-ray observations 
on electron temperature profiles and with some theoretical and numerical 
predictions about mean temperature profiles. 



Subject headings: galaxies: clustering — galaxies: intergalactic medium - 
galaxies: X-ray — hydrodynamics 



1. INTRODUCTION 



Clusters of galaxies (CG) are the largest virialized objects in the universe, which contain 
collisionless particles, galaxies and dark matter, and a diffuse gas component. The gas 
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component is called intracluster medium (ICM). The ICM is the plasma with temperature 
of about 10 8 K, thus emitting X-ray mainly through the thermal bremsstrahlung of the 
electrons (Sarazin 1988). 

Revealing temperature profiles of ICM is an important problem. In theoretical work 
concerning CG, isothermality of ICM is often assumed. In particular, the isothermal (3 
model (Cavaliere k Fusco-Femiano 1976) is adopted in the conventional mass determination 
through X-ray observations, in the estimation of the Hubble constant (Hq) through 
Sunyaev-Zel'dovich (SZ) effect (Sunyaev k Zel'dovich 1972; Birkinshaw, Hughes, k Arnaud 
1991), etc. However, deviation from a flat temperature profile becomes important at lower 
densities in the mass estimation using X-ray data (Evrard, Metzler, k Navarro 1996; 
Schindler 1996). Also in the H estimation through SZ effect the non-isothermality is 
one important source of errors (Inagaki, Suginohara, k Suto 1995; Yoshikawa, Itoh, k 
Suto 1998). Therefore, accurate temperature profiles of ICM are required to improve such 
methods. 

On the other hand, temperature maps of ICM provide us with useful information 
about CG. In merging clusters, characteristic temperature structures are expected to 
occur through shock heating and adiabatic compression. Some numerical simulations are 
especially focused on this problem (Schindler k Miiller 1993; Ishizaka k Mineshige 1996; 
Roettiger, Loken, k Burns 1997; Ishizaka 1998; Roettiger, Stone, k Mushotzky 1998). In 
particular, off-center collisions are investigated by Ricker (1998). Comparing these results 
with the X-ray observational data , we can guess in what phase the merging clusters are 
(Fujita et al. 1996; Honda et al. 1996; Churazov et al. 1998; Donnelly et al. 1998; Davis k 
White 1998). 

Recent X-ray observations with ASCA and ROSAT reveal the radial electron 
temperature distribution of ICM (Markevitch et al. 1996; Markevitch 1996; Markevitch, 
Sarazin, k Irwin 1996; Markevitch et al. 1997). In some clusters the observed electron 
temperature gradients correspond to the polytropic index, j(P oc p 7 ), of 7 = 1.5 or 
even more (Markevitch 1996). On the other hand, 7 ~ 1.2 is expected for plasma mean 
temperature distribution derived by the self similar solution of Bertschinger (1985) and 
some numerical simulations (Evrard 1990; Katz k White 1993; Navarro, Frenk, k White 
1995; Eke, Navarro k Frenk 1997; Takizawa k Mineshige 1998; Bryan k Norman 1998). 
There is a discrepancy between the observed values and the theoretical expectations. 

To explain the discrepancy between the observed electron temperature profiles and 
the theoretical suggestions on the mean temperature profiles consistently, we construct a 
model of ICM incorporating properly the relaxation process between the ions and electrons. 
Fox k Loeb (1997) was the first to investigate the two-temperature nature of ICM. They 
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construct the electron and ion temperature profiles by combining the self-similar solution 
by Bertschinger (1985) with the analytical evolutionary model of adiabatic two-temperature 
plasma which is originally found by Shafranov (1957). Chieze, Alimi & Teyssier (1998) 
carried out the 3-D hydrodynamical simulations of two-temperature ICM. However, CGs 
only in the Einstein de Sitter universe were considered in the previous work concerning 
the two-temperature model of ICM. Since in a higher density universe CG form at more 
recent epoch (Richstone, Loeb, & Turner 1992), the previous work is restricted to the 
case in which a temperature difference is expected to be the largest among reasonable 
cosmological models. Therefore, it is necessary to study a two-temperature model of ICM 
in other cosmological models in order to confirm whether the temperature difference is 
really significant or not. This problem is also related to the dependence of temperature 
profiles on the cosmological parameters, which is discussed in some numerical simulations 
(Evrard et al. 1996; Eke et al. 1997). Furthermore, quantification of resultant temperature 
profiles by using polytropic indices is not fully discussed in the previous work, which is very 
important to compare the models with X-ray observations properly. Ettori & Fabian (1998) 
studied two-temperature ICM using simple analytic models and applied their results to the 
CG A2163. However they neglected the dynamical properties of ICM and assumed that all 
the ICM in CG is heated at the same time. This assumption is problematic to quantify the 
temperature profiles because the timescale of shock propagation in CG is comparable to the 
dynamical timescale and not much shorter than the age of the universe. 

For these purpose, we apply the method of Fox & Loeb (1997) to the results of 
numerical simulations instead of the self-similar solution. Thus we can investigate 
two-temperature ICM quantitatively in cosmological models other than the Einstein de 
Sitter model. Studying the dependence of electron temperature profiles on cosmological 
models is another important purpose. 

The rest of this paper is organized as follows. In §|2| we estimate the relevant timescales 
and spatial scales of ICM. In §|3| we briefly review the adiabatic model of two-temperature 
plasma proposed by Fox & Loeb (1997). In §|| we described the method to calculate the 
ion and electron temperature profiles of ICM by combining a one-temperature ICM model 
with the adiabatic model by Fox & Loeb (1997). In §[5] we describe the adopted numerical 
methods and initial conditions for simulations of CG. In §|6| we present the results. In we 
summarize the results and discuss their implications. 
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2. ORDER ESTIMATION 

We consider fully ionized plasma which consists of electrons and ions. Only coulomb 
coupling is considered as the relaxation process. Then two-body relaxation timescale of a 
x-particle, whose density is n x and temperature is T x , is (Spitzer 1962) 

m^(3kT x )^ 

xx 5.7l7m x e 4 Z x MnA' 1 } 

where, m x is the particle mass, Z x is the particle charge number, e is the electron charge, k 
is the Boltzmann constant, and In A is the Coulomb logarithm and approximated to be 

lnA~37.8 + ln(4M--ln( "^—^ (2) 

U0 8 K/ 2 V10- 3 cm- 3 / v ; 

for T x > 4 x 10 5 K. Therefore, t„ is larger than t ee by a factor of the order of [m-Jmc) 1 ^ 2 . 

On the other hand, the equilibrium timescale between ions and electrons is, 

3m c mj / kT c kTj \ 3 / 2 

ei ~ 8(27r) 1 /% i Z i 2 e 4 InA Vm c m-J ' 1 ] 

Therefore, t ei is grater than tn by a factor of (m i /m c ) 1 '' 2 . In addition, this timescale can be 
comparable to or longer than the Hubble time in the outer region of CGs, 

y V 40 J vio- 3 cm- 3 / vio 8 k; v ; 

Suppose that ICM is heated through the shock in accretion flow (Cavaliere, Menci, 
& Tozzi 1997). The shock primarily heats ions because the kinetic energy of a particle 
is proportional to the particle mass. In the post shock region the ions reach thermal 
equilibrium on a timescale of tn after they are heated through the shock. Within this time 
ion temperature is significantly higher than electron one. Eventually thermal energy is 
transported from the ions to the electrons through the Coulomb collisions between the ions 
and electrons and 7] becomes comparable to be T e on timescale of t ei . Note that t ci > t ce , tn. 

Under such circumstance, the radial length, r tt , over which the electron temperature 
is significantly lower than the ion temperature can be estimated as follows. Denote the 
propagation speed of the shock front by v shock- Then, we have, 

rtt - tehock- (5) 

Using the strong shock approximation and neglecting the post-shock gas velocity compared 
with v shock, we find 

fshock - gfinfall, (6) 
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where tWaii is the infalling velocity of the gas, which is related to the post-shock gas 
temperature, T. The kinetic energy of infalling gas is nearly equal to the thermal energy of 
the post shock gas; namely 



^P^nfall ~ k BT. (7) 



Therefore, using the equations (f|), (||), (|6]), and (|7|), we derive 

1 V 40 J V10- 3 cm- 3 / V10 8 K7 v ; 

Importantly this is smaller by about one order of magnitude than the spatial scale of CGs. 
Note that the higher temperature is, the wider becomes the region where the temperature 
difference between ions and electrons is significant. 

Another important timescale related to ICM is the radiative cooling timescale, t c . 
Since thermal bremsstrahlung is the dominant cooling process in ICM, we have 

t c = 8.5 x 10 10 yrf ?^Vf Y -^-Y^ (9) 

y V10- 3 cm- 3 / V10 8 K/ v ; 

From equations (|J) and @ we find, 

^ = 2.3xl0-^r(4-Y (10) 

t c v 40 J vio 8 k; v ; 

Hence t c is always longer than t e i in the typical ICM, and thus we can safely neglect cooling 
effects as far as we are concerned with the overall cluster structure. 

On the other hand, the age of CG, t age , is of the order of 10 9 or 10 10 yr. Therefore, as 
long as heating from galaxies can be neglected, we can divide intra-cluster space into three 
regions according to the magnitudes of these three timescales, t ei , t c , and t age as follows: 



1. The central higher density region where t e i < < ^ age - Radiative cooling is important 
but ICM can be regarded as one-temperature fluid. This situation corresponds to the 
so-cold 'cooling flow'. 

2. The middle region where t e i < ^a g0 < t c - We can regard ICM as adiabatic 
one-temperature fluid. 

3. The outer lower density region where t age < t ei < t c . Radiative cooling can be 
negligible and the electron temperature is considerably lower than the ion temperature. 
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3. EVOLUTION OF ADIABATIC TWO-TEMPERATURE PLASMA 

The formulation here is based on Fox & Loeb(1997). We neglect thermal 
conduction, which is the case if tangled magnetic fields suppress conduction. The 
Laglangean time evolution of the electron temperature, T e , and the mean temperature, 
T = (n e T e + mTi)/(n e + ni), in the adiabatic fluid element is, 

at t e i n at 

dT , . T dn .„„. 

* = ^- 1] nTV (12) 

where TJ is the ion temperature, n is the gas density, and 7 = 5/3, is the ratio of specific 
heat. Introducing the temperatures normalized by T, T e = (T e /T) and T ; = (T;/T), we find, 

dT e T\ T e . . 

~ 3/2 

Note that £ e i is proportional to T e In A, since in the post-shock region the gas behaves 
adiabatically (T oc n 2//3 ). Thus t 2s = t ei (t)T c (t)~ 3 / 2 is constant in time, if we neglect the 
small change due to the Coulomb logarithm. Now equation ([□!) becomes 

dT e 1 fn\-\-n e \, ~. --3/2 , <N 

/ 1 • W l - T e )T e . (14) 



dt t 2s \ n 



Since ICM is almost perfectly ionized, the ratio of (nj + n c )/ni can be regarded as constant. 
Thus equation (JTJ) can be integrated analytically. If we assume the pre-shock T e equal to 
zero, the solution for the fluid element which has passed the shock front at t — t 8 is, 



t — t s — t2s 



rii + n c 



1 - \ T t 



(15) 



Thus we can obtain the Laglangean time evolution of T e of the fluid element after the 
passage through the shock by solving equation (|15|). 



4. NUMERICAL METHOD 

Using equation ([15]) we can construct the temperature profiles at t — t , T e (r,t ) and 
Ti(r,t ), as follows. First, we give the velocity field of the gas, v(r, t), the shock radius, 
fshock (t), and the equilibrium timescale of the fluid element at the passage through the 
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shock front, t 2s [r shockit)}, using some analytical models or numerical simulations presented 
later. Then we can calculate the Laglangean path of the fluid element which passed the 
shock surface at t — t s , R(t; t B ), by solving the ordinary differential equation, 

dR 

— = v(R(t;t s ),t), (16) 

with the appropriate initial condition, 

R(t s ;t B ) = r shock (t s ). (17) 

Integrating differential equation (|T6| ) for various t s , we obtain r = R(t ;t s ), the radius 
which the fluid element that passed the shock at t — t s resides at t — to, as a function of t s . 
Thus we can regard to — t s as a function of the radius, 

t -t a = f{r;t ). (18) 

Using the equations (p~5|) and (0), and the model of ^[^shock^s)], we can solve T e at t = t 
as a function of r. Finally, we obtain T e (r) and T[(r) using T e (r) and the model of T(r). 



5. THE SIMULATIONS 

5.1. Numerical Method 

To give f(r, t), r s h oc k(£), etc, we have performed numerical simulations of a spherically 
symmetric CG. For dark matter (DM), we use the shell model (Henon 1964). We set the 
number of shells, N, equal to 5000. As for gas on the other hand, we use 1-dimensional, 
spherically symmetric, total variation diminishing (TVD) code with minmod limiter (Hirsch 
1990). Note that TVD code is one of the most powerful tools to treat shocks. One mesh 
spacing corresponds to Ar = 0.005/(1 + z) Mpc. We assume that the gas is ideal, with 
7 = 5/3. As to the boundary conditions, the inner edge is assumed to be a perfectly 
reflecting point. The outer edge is assumed to be a perfectly transmitting surface. The 
basic equations and the numerical method used here are fully described in §2 of Takizawa 
& Mineshige (1998). 



5.2. Models and Initial Conditions 



In this paper, all of the calculations are carried out from z ini = 10 to the present 
time (z = 0). The cosmological models which we adopt are (ri , A ) = (1.0,0.0) (ED), 
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(fi ,A ) = (0.2,0.0) (OP), and (fi ,A ) = (0.2,0.8) (FL). The Hubble constant, H is set 
to be Hq = 100kms _1 Mpc~ 1 in the transformations of length and time coordinates. We set 
fib = 0.012/i -2 in all models taken from the nucleosynthesis determination (Walker et al. 
1991). Note that fi DM = fio - fib- 

We make initial density profiles in the same manner as Peebles (1982). At first we 
prepare N concentric shells with a constant density equal to fioM at z ini = 10. Then, 
a density fluctuation is introduced by perturbing the radius and velocity of each shell 
following 



where rf is the unperturbed coordinate, S(r) represents the mean density fluctuation inside 
r, which is derived from the density fluctuation field (specified below) within r, and H(zi n i) 
is the Hubble constant at z — z^. Here, we used the Zel'dovich approximation (Zel'dovich 
1970) and the approximation of dlogDi/dloga ~ fi 06 , where D\ is the linear perturbation 
growth rate for the growing mode and a is the scale factor (see Suto 1993). 

The initial conditions of the density fluctuation field, S(r), are generated by applying 
the Hoffman- Ribak method (Hoffman & Ribak 1991; van de Weygaert & Bertschinger 
1996) to spherical systems (see Appendix). We constrained in such a way that there exist 
initially density enhancements on 1 Mpc scale whose amplitudes correspond to 3<j level in 
the CDM power spectrum. The normalization is og = 1 in each cosmological model, which 
was obtained from the observation of the nearby galaxy distribution (see Suto 1993). 

The initial conditions of gas are set as follows. At first, gas density was everywhere 
taken to be the mean baryon density of the universe at z — 10, and the temperature of the 
gas (T gaSj i) was constant everywhere; T gaS)i = 10 7 K in all models. We then add the adiabatic 
fluctuation in such a way that the ratio of the DM density and the gas density remains the 
same. Note that after the perturbation is added, the gas temperature distribution becomes 
nonuniform, accordingly. Moreover, note that temperature of the infalling gas at z ~ 1 is 
sufficiently lower than virial temperature, since the gas expands adiabatically following the 
cosmological expansion. 




(20) 



(19) 
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6. RESULTS 

The overall evolution of the simulated CG is essentially the same as that presented 
by Takizawa & Mineshige (1998). Thus we concentrate on showing the results of the 
temperature profiles, mass estimation, and so on. 



6.1. Temperature Profiles 

Figure [I] shows radial distribution the T e (r) (solid line), 71(r) (dotted line), and T(r) 
(short dashed line) of Model ED at z = 0. At r < 0.6Mpc T e is very close to 7], while, at 
r > 0.6Mpc the electron temperature is considerably lower than the ion temperature and 
the discrepancy gets increased outward. At r ~ IMpc, the electron temperature is only a 
half of the mean temperature. Small scale fluctuations in the temperature profiles are due 
to the sound wave propagation in ICM (Takizawa & Mineshige 1988). 

In real X-ray observations what is actually obtained is an emissivity-weighted, 
line-of-sight projected electron temperature map, which is displayed in Figure (^) for model 
ED. In this figure, we assumed the spatial resolution to be 0.25 Mpc, which corresponds to 
about 6' for an object located at z = 0.05, and the error of the temperature measurement 
to be 20 %. Also in this map we can clearly see that the electron temperature decreases 
outward. 

To estimate the temperature gradients both for T e and T, we measure the polytropic 
indexes, 7 P , in the usual way. We fit the density profile by the /3-model and the temperature 
profile by the polytropic model as follows, 



n(r) = n Q 



1+ 



r N 2-| -3/3/2 



(21) 



T(r) oc n(r) 7 ^ 1 (22) 

where no, r c , and j3 is the fitting parameters of the /3-model. We fit the resultant density 
profile only inside the shock front. The data are fitted by chi-square fitting. We assume 
that the variance in any quantity is proportional to its square (Sf/f = const.) because our 
main purpose is not to simulate observations with a specific instrument but to obtain the 
intrinsic profiles of our calculated results. The results of the fitting are summarized in Table 
[I] and the corresponding polytropic indices are listed in Table |2|. From Table [1| we find 
that the f3 values are somewhat bigger than those typically observed. On the other hand, 
our results are well coincident with the self-similar solution by Bertschinger (1985). Thus 
the bigger (3 values are probably due to the assumption of spherical symmetry and neglect 
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of angular momentum of gas. If the initial gas temperature is higher, which corresponds 
to the case that reheating of the ICM, e.g., by protogalaxies, is substantial, f3 values can 
be smaller (Metzler & Evrard 1994; Takizawa & Mineshige 1988). From Table |2| we find 
that the polytropic indices of the electron temperature profiles are 7 P ~ 1.5, which are 
systematically larger than those of the mean temperature profiles, 7 P ~ 1.3. 

The specific entropy profiles derived from the electron temperature [S e oc m(T e /p 7-1 )], 
and from the mean temperature [S oc ln(T/p 7-1 )], are shown in Figure [5] by the solid and 
dotted lines, respectively. Entropy is normalized to be zero at the inner boundary. The 
latter S rises outward, which is characteristic of the ICM heated through the shock (Evrard 
1990; Takizawa & Mineshige 1998). The former S e , in contrast, rises outward as the case 
only in the inner region (r < 0.6Mpc), stays nearly constant in the middle region, and falls 
outward in the outer region. 

6.2. Dependence of Fitting Parameters on Size of Fitting 

We fit the calculated density profile inside the shock front by equation Q2"I|) and listed 
the fitting results in Tables [l] and |[ However, it is possible that the fitting results can be 
influenced by the position of the outer edge of the region used for the fitting. In real X-ray 
observations the outer edge of the X-ray emitting region is perhaps inside the shock front. 
To assess this effect we fit the density and temperature profiles inside various radii for 
Model ED and list the dependence of the fitting results related to the density profile (no, 
r c , and 0) on the outer-edge radius (r out ) in Figure £|. It is found that these parameters are 
insensitive to the outer radius as long as r out > 0.8 Mpc. In density profiles, therefore, we 
can safely neglect the influence of the outer edge. This fact is actually expected because of 
the self-similar nature of the gas density profile (Bertschinger 1985; Takizawa & Mineshige 
1998). 

On the other hand, the dependence of 7 P on r out is rather different. Figure shows 
that 7 P derived from T e (by closed squares) monotonically increase outward whereas 7 P 
derived from T (by open square) does not exhibit systematic changes. When we calculate a 
polytropic index from X-ray observational data, thus, the resultant value could be subjected 
to large errors arising from the finite detection limit, background noise, and so on. 

The results of model OP and FL are essentially the same as those of model ED. 
Although there is a simple self-similar solution only for model ED, self-similar nature is also 
expected both in OP and FL (Takizawa & Mineshige 1988). 
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6.3. Mass Estimation 

Since the observed electron temperature significantly deviates from the mean 
temperature that determines the dynamics of the system, the total mass of CG is probably 
underestimated if hydrostatic equilibrium is calculated based on the electron temperature 
map. When ICM is assumed to be isothermal we usually use the emissivity-weighted mean 
temperature, which is more like the temperature in the central high-density region. Thus, 
the underestimation of the mass is practically negligible. When the electron temperature 
profile is used, conversely, the mass can be seriously underestimated. Hence the mass 
derived from the assumption of hydrostatic equilibrium, M hydro (r), is, 

M hydlo [r } - ^ G r [x + (r/rc)2]1+3/3(7p _ 1)/2 - W 

Figure |6| depicts the ratio of Mbydro to the actual mass as a function of radius. The solid 
line represents the ratio calculated based on T e (r) and the dotted line represents that based 
on T(r). The mass derived from T e is underestimated by almost 50 % because of the lower 
electron temperature. Note that the mass derived from T is also underestimated by about 
10 % due to the bulk motion of the gas. 

In the central region Mh y dro based on T e and T are both larger than the actual mass. 
This is due to the fitting errors of temperature profiles. Since the core radii are different 
between the density and temperature profiles, polytropic model cannot describe the our 
results well in this region. 



7. SUMMARY AND DISCUSSION 

We constructed the models of ICM, incorporating the relaxation process between the 
ions and electrons. From the simple order estimation, we find that the electron temperature 
is well below the ion temperature in the outer region of CG and that such a lower T e region 
spreads over a Mpc scale in typical CG. In addition, the hotter CG is, the wider becomes 
the two-temperature region. Comparing three relevant timescales in ICM (the age of CG, 
radiative cooling, and equilibrium timescales between ions and electrons), we can divide 
ICM into three regions; from the center outward, the cooling dominant, one-temperature 
region, adiabatic, one-temperature region, and adiabatic, two-temperature region. 

We calculate the temperature profiles of two-temperature ICM combining the spherically 
symmetric, N-body and hydrodynamic simulations for three different cosmological models 
with the adiabatic two-temperature plasma model by Fox & Loeb (1997). While the 
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polytropic indices of the mean temperature profiles are ~ 1.3, those of the electron 
temperature profiles are ~ 1.5. As a consequence, the specific entropy profiles derived from 
the electron temperature are rather flat. 

We examine the dependence of the fitting parameters on the radius of the outer edge 
of the region used for the fitting when the density profile is fitted by the /3-model and the 
temperature profile by the polytropic model. The fitting results of the density profile and 
mean temperature profile is insensitive to the outer edge. On the other hand, the polytropic 
index derived from the electron temperature rises as the outer-edge radius increases. 

The total mass of CG is underestimated about 50 % when we use the electron 
temperature profile though the underestimation is negligible when we assume that ICM is 
isothermal and adopt emission-weighted mean temperature as the temperature of ICM. 

We confirm that temperature difference between ions and electrons in ICM is 
substantial in cosmological models other than the Einstein de Sitter model. The polytropic 
indices of electron temperature profiles are insensitive to the cosmological models in the 
range of our calculations. Note that the baryon density of these models is set to be the 
same. If the baryon fraction is set to be constant in each model, the result will probably 
changes because the equilibrium timescale is sensitive to the baryon density. 

In general our resultant electron-temperature profiles tend to have steeper gradients 
than those of Ettori & Fabian (1998). They assumed that all the ICM in CG is heated 
at the same time and that ICM is in hydrostatic equilibrium. On the other hand, in our 
calculations the ICM in the outer region is heated more recently than that in the inner 
region because the shocks propagate at finite speed. Furthermore, gradual radial infall 
persists inside the shock fronts and makes the two-temperature region slightly compress 
inward. (Takizawa & Mineshige 1998). Thus the two-temperature region of our results 
becomes wider than those of Ettori & Fabian (1988). 

Although only spherical accretion is considered in this paper, there may also arise 
asymmetric merging of comparable clumps in reality. It is possible that shocks occurring 
in merging events will also generate the temperature difference between ions and electrons. 
In this case it is believed that a bow shock with an arc shape is formed just between 
the centers of two substructures (Schindler & Miiller 1993; Ishizaka & Mineshige 1996; 
Roettiger et al. 1997; Ishizaka 1998). In the post-shock region the energy is transported 
from ions to electrons but that timescale changes along the shock front since the timescale 
is sensitive to the density of ICM. Therefore, if we consider the temperature difference 
between ions and electrons, the location of the observed hot gas region is probably shifted 
and the shape is more deformed in comparison with the results of the former simulations. 
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We will investigate this issue as a future work. 

We consider only the classical coulomb coupling as the relaxation process between 
ions and electrons. It is possible, however, in ICM more efficient relaxation processes may 
be effective (McKee & Cowie 1977; Pistinner, Levinson, & Eichler 1996). In this case the 
equilibrium timescale can be shorter than the value given by equation (f|). Therefore, the 
temperature difference between ions and electrons can be less and polytropic indices can 
be smaller than our results. If magnetic field exits in ICM, it is possible that electrons are 
also significantly heated in shocks by MHD instabilities. Also in this case the temperature 
difference between ions and electrons can be less. 

We neglect the heating process from the galaxies to ICM. If thermalized hot gas is 
injected to ICM from the galaxies, temperature difference could be less. Furthermore, in 
this case the preheating of the ICM influences the density profile and mean temperature 
profile (Metzler & Evrard 1994; Navarro et al. 1995; Takizawa & Mineshige 1998). Since 
heavy elements like iron are detected in ICM at redshift up to z ~ 1 (Hattori et al. 1997), 
this effect should be considered to construct a more realistic model. In such a study the 
adiabatic model like Fox & Loeb (1997) cannot be adopted. Thus non-adiabatic models or 
fully two- fluid simulations like Chieze et al. (1998) are required. 

The author would like to thank S. Mineshige for valuable comments. This work is 
supported by a Grant-in-Aid from the Ministry of Education, Science, Sports and Culture 
of Japan (6179) and Research Fellowships of the Japan Society for the Promotion of Science 
for Young Scientists. 

A. THE HOFFMAN-RIB AK METHOD FOR SPHERICAL SYSTEMS 

We consider a random homogeneous and isotropic Gaussian field / with zero mean 
which is defined by its power spectrum P(k). When / is subjected to linear constraints, 
g = Cf, then the constrained Gaussian field / is realized as follows (Hoffman & Ribak 



where / is an unconstrained random Gaussian field whose power spectrum is P(k), M is 
the two-point correlation matrix obtained from P(k), and Q = CMC^ . 

In general three dimensional case, / is described by its Fourier components /k, 



1991) 



f = f + MC^Q-\g-Cf) 



(Al) 




(A2) 
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However, to generating the initial conditions which can be used for our spherical 
symmetric simulations, we should construct / as a function of only r instead of r. Since / 
is isotropic f(r) is obtained as follows. In spherical systems it is convenient to expand /(r) 
in spherical harmonics times spherical Bessel functions ji (Binney & Quinn 1991). A plain 
wave is expanded in these functions as follows, 

+00 +1 

exp(*kr)=47r£ £ i l ji(kr)Yr(O k ,<f> h )Yr(e r ,<i> r ). (A3) 

1=0 m=-l 

Integrating this for the tangential components we find, 

J exp(ikr)dk = Ait J j (kr)k 2 dk. (A4) 

We can assume that /k = fk because of isotropy of /. Thus we obtain, 
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Model 


n (10 3 cm 3 ) 


r c (Mpc) 


P 


ED 


1.16 


0.067 


0.793 


OP 


3.64 


0.104 


0.815 


FL 


1.45 


0.155 


0.903 



Table 1: Density profiles of gas at z — 0. 
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Model 


7 P of T c 


7 P of T 


ED 


1.55 


1.35 


OP 


1.44 


1.34 


FL 


1.51 


1.37 



Table 2: Polytropic indices of T c and T at z = 0. 
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Fig. 1. — The profiles of the T e (solid line), T ; (dotted line), and T (short dashed line) of 
Model ED at z = 0. At r < 0.6 Mpc T e ~ T h whereas at r > 0.6 Mpc T c is considerably 
lower than 7] and the discrepancy gets enhanced outward. At r ~ 1 Mpc, T c ~ T/2. Small 
scale fluctuations in the temperature profiles are due to the sound wave propagation in ICM. 

Fig. 2. — Emissivity weighted, line-of-sight projected electron temperature profile of Model 
ED at z — 0. Here we assumed the spatial resolution to be 0.25 Mpc, which corresponds to 
about 6' for an object located at z — 0.05, and the error of the temperature measurement 
to be 20 %. Clearly the electron temperature decreases outward. 

Fig. 3. — The specific entropy profiles derived from the electron temperature [S e oc 
ln(T e /p 7-1 )], by the solid line and from the mean temperature [S oc ln(T/p 7-1 )], by the 
dotted line, respectively. Entropy is normalized to be zero at the inner boundary. 

Fig. 4. — The dependence of the parameters related to the density profile (n , r c , and (3) 
on r out , the radius of the outer edge of the region used for fitting. It is found that these 
parameters are insensitive to the outer edge radius as long as r out > 0.8. 

Fig. 5. — The dependence of 7 P on r out . We find that 7 P obtained from T e (closed squares) 
increases as r out increases, whereas 7 P obtained from T (open squares) is insensitive to r out . 

Fig. 6. — The ratio of Mh y dro to the actual mass. The solid line is the result based on T e {r) 
and the dotted line is that based on T(r). The mass derived from T e is underestimated by 
almost 50 % because of the lower electron temperature. Note that the mass derived from T 
is also underestimated by about 10 % due to the bulk motion of the gas 



3 — 







I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



2 — 



1 — 



I E- 



i i i i i i i i i i i i i i i i i i i 



— EE 1 i 

i i i I i i i I i 



.2 .4 



.6 .8 
r [Mpc] 



1 1.2 1.4 





9 




1.5 






'e 




o 




1 


1 


o 




o 
C 


.5 




.2 

.15 



o 

Q. 
5, -1 



i r 



i i r 



i i r 



i r 



□ 



□ 



□ 



□ 



□ 



(a) 

□ □ 



J L 



J I L 



J I L 



J L 



.8 



Tout [Mpc] 



.1 



1.2 



i r 



i i r 



i i r 



i r 



(b) 



□ 



□ 



□ 



□ 



□ 



□ 



.05 




J L 



J I L 



J I L 



□ 



J L 



.6 



.8 



r out [Mpc] 



,1 



1.2 





I 


I 


I 


I 


I 


I 


I 

(c) 


I 




□ 


□ 
















□ 


□ 


□ 


□ 


□ 






I 


I 


I 


I 


I 


I 


I 


I 



^ 1 

.5 




.8 



r ou t [Mpc] 



.1 



1.2 



1.6 



1.5 



£ 1.4 



1.3 



1.2 





I 


I 


I 


I 


1 


1 


1 


1 


- 














■ 


— 















1 















1 


_ 
■ 














_ 
■ 












— 


■ 


■ 


□ 


□ 








— 






n 

I I 












□ 


































□ 


□ 


□ 




- 










— 




I 


I 


I 


1 


1 


1 


1 


1 



.6 .8 1 1.2 

Tout [Mpc] 



